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Repeated psychosocial stress in early-life has significant impact on both behavior and neural 
function which, together, increase vulnerability to depression. However, neural mecha- 
nisms related to repeated stress remain unclear. We hypothesize that early-life stress may 
result in a reduced capacity for cognitive control in response to a repeated stressor, partic- 
ularly in individuals who developed maladaptive emotional processing strategies, namely 
trait rumination. Individuals who encountered early-life stress but have adaptive emotional 
processing, namely trait mindfulness, may demonstrate an opposite pattern. Using a men- 
tal arithmetic task to induce mild stress and a mindful breathing task to induce a mindful 
state, we tested this hypothesis by examining blood perfusion changes over time in healthy 
young men. We found that subjects with early-life stress, particularly emotional abuse, 
failed to sustain neural activation in the orbitofrontal and ventromedial prefrontal cortex 
(vmPFC) over time. Given that the vmPFC is known to regulate amygdala activity dur- 
ing emotional processing, we subsequently compared the perfusion in the vmPFC and 
the amygdala in depression-vulnerable (having early-life stress and high in rumination) 
and resilient (having early-life stress and high in mindfulness) subjects. We found that 
depression-vulnerable subjects had increased amygdala perfusion and reduced vmPFC 
perfusion during the later runs than that during the earlier stressful task runs. In contrast, 
depression-resilient individuals showed the reverse pattern. Our results indicate that the 
vmPFC of depression-vulnerable subjects may have a limited capacity to inhibit amygdala 
activation to repeated stress over time, whereas the vmPFC in resilient individuals may 
adapt to stress quickly. This pilot study warrants future investigation to clarify the stress- 
related neural activity pattern dynamically to identify depression vulnerability at an individual 
level. 



Keywords: early-life stress, repeated stress, depression vulnerability, depression resilience, fMRI, ventromedial 
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INTRODUCTION 

Stress is a significant risk factor for depression and anxiety. Chronic 
stress can produce significant detrimental effects psychologically 
and physiologically (McEwen and Stellar, 1993). In response to 
chronic stress, some individuals develop neural habituation and 
adapt to the stress, but others may become sensitized to the stres- 
sor and may have prolonged or amplified neural responses over 
time (McEwen and Stellar, 1993). Studying neural response to 
a repeated stressor over time enables us to characterize neural 
signatures to identify depression vulnerability. The pattern of indi- 
vidual differences in response to a repeated stressor over time 
may become an important feature in differentiating vulnerabil- 
ity to depression and anxiety from resilience at an individual level. 



However, neuroimaging studies on repeated stress are rare, and the 
impact of repeated stress on neural circuits related to depression 
vulnerability remains unclear. 

McEwen and Stellar (1993 ) proposed four types of maladaptive 
responses to repeated stress: (1) lack of habituation, (2) pro- 
longed response, (3) inability to reduce reactivity, (4) inability 
to respond. For the first three types of responses, stronger neural, 
or physiological responses may be found in response to a repeated 
stressor compared to a single-time stressor. Studies in peripheral 
reactivity such as blood pressure have demonstrated this reac- 
tivity pattern (Schneider et al., 2003). However, brain regions 
that reflect this response pattern to repeated stress are still under 
investigation. 
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By contrast, a number of brain regions such as the ventro- 
medial prefrontal cortex (vmPFC), ventrolateral prefrontal cor- 
tex (vlPFC), amygdala, and hippocampus can reveal maladaptive 
response patterns in individuals vulnerable to stress in response 
to repeatedly presented stressors (McEwen, 1999; Taylor et al., 
2004, 2011). Dysfunction of these regions while processing emo- 
tional events are frequently reported in individuals with depression 
(Siegle et al, 2007; Fitzgerald et al, 2008; Wang et al, 2008), anx- 
iety, post-traumatic stress disorder (PTSD) (Shin et al., 2004a,b), 
and in individuals vulnerable to depression such as those with 
an early- life stress history (Matsumoto et al, 2009; Pechtel and 
Pizzagalli, 2011). Early-life stress history refers to exposure to 
emotional or physical abuse, such as neglect and harsh and 
chaotic parenting in childhood (Bernstein et al., 2003). Both the 
vmPFC and vlPFC have direct projections to the amygdala (Sah 
et al., 2003; Holland and Gallagher, 2004) and have been found 
to exert a top-down, inhibitory effect on the amygdala. Typi- 
cally, while regulating negative affect, healthy individuals have 
demonstrated higher vmPFC (Urry et al, 2006) or vlPFC (Tay- 
lor et al., 2006) activation and lower amygdala activation. A study 
on early-life stress (Taylor et al, 2006) found that in response 
to threatening cues, healthy controls showed a negative correla- 
tion between activation of the right vlPFC and amygdala activ- 
ity, whereas participants with an early-life stress history showed 
a positive correlation between the two regions. These results 
suggest a possible failure to recruit the prefrontal cortex effec- 
tively for regulating emotional responses to threatening cues, and 
thus no reduction in amygdala activity was found in individuals 
with early-life stress. Therefore, it is important to further clar- 
ify the neural response pattern between the prefrontal cortex and 
the amygdala to repeated stressors in individuals with early-life 
stress. 

We hypothesized that individuals with early-life stress would 
show an activation pattern in response to a repeated stressor 
that reflects long-term stress exposure, characterized by sustained 
amygdala hyperactivation and attenuated activation in the vmPFC 
or vlPFC over time due to failed recruitment of sustained activ- 
ity from the prefrontal cortex to regulate amygdala response 
continuously. 

While history of early-life stress undoubtedly is a risk factor 
for future depression (Wals and Verhulst, 2005), not everyone 
with early-life stress eventually develops depression (Charney and 
Manji, 2004). Individual differences in the ability to regulate 
stress and emotional reactions to stressful situations can mod- 
ulate brain activity and peripheral physiology, with implications 
for vulnerability vs. resilience to depression. For example, indi- 
viduals who tend to ruminate about chronic stress and negative 
emotions develop depression more easily in response to acute 
life stress (Nolen-Hoeksema, 1991). Rumination is defined as 
"repetitively focusing on the fact that one is depressed; on one's 
symptoms of depression; and on the causes, meanings, and conse- 
quences of depressive symptoms" (Nolen-Hoeksema, 1991). It has 
been consistently found that rumination increases negative mood 
(Ward et al, 2003) and predicts the onset and relapse of depres- 
sion (Nolen-Hoeksema, 1991; Nolen-Hoeksema et al, 1993; lust 
and Alloy, 1997). Therefore, rumination is a key cognitive trait 
underlying vulnerability to depression. 



Some individuals with early-life stress may develop traits that 
enable them to be resilient (Marks et al., 2010; Binder et al., 2011; 
Carli et al, 2011). One such trait is mindfulness. Trait mind- 
fulness, also known as dispositional mindfulness, refers to the 
self-regulation of attention as well as an orientation of open- 
ness, curiosity, and acceptance to all experiences (Bishop et al, 
2004). Individuals who are high in trait mindfulness in daily life 
demonstrate better psychological health and lower incidence of 
depression (Keng et al., 2011). Therefore, it is conceivable that 
those who have developed trait mindfulness in spite of having 
an early-life stress history have a different neural activity pattern 
in response to a repeated stressor compared to individuals with 
early-life stress who have developed trait rumination. 

We hypothesized that individuals with early-life stress, will 
show increasingly or consistently high activation in the amyg- 
dala and low or quickly reduced cognitive regulatory responses of 
the vmPFC following a repeatedly presented mild stressor during 
later relative to early trials. We predicted this activation pattern to 
repeated stress over time will be found particularly in individuals 
with early-life stress and high trait rumination, whereas an oppo- 
site pattern will be found in individuals with early-life stress and 
high trait mindfulness. 

MATERIALS AND METHODS 
PARTICIPANTS 

To avoid the confounding influence of menstrual cycle on stress 
sensitivity (Ossewaarde et al., 2010), only male participants 
[n = 19, mean (SD) age = 27.05 (7.21)] were recruited in the study. 
Participants were recruited from the subject registry at the Duke- 
UNC Brain Imaging and Analysis Center. The exclusion criteria 
were: (1) MRI contradictions, (2) current or history of neuro- 
logical and psychiatric disorders, (3) drug abuse, and (4) current 
medication use. The study was approved by the Duke University 
Health System Institutional Review Board. All participants have 
signed written consent. 

PROCEDURES 

A stress induction task and a mindful breathing task were admin- 
istered over two different days separated by 7-10 days. The 
order of the stress and mindfulness tasks was counterbalanced 
among the participants. Each day was composed of a pre-scan 
stress/mindfulness task training session and an imaging scan ses- 
sion. Detailed procedures can be found in our previous report 
(Paul et al., 2013). A number of questionnaires were adminis- 
tered during pre-scan session to quantify affective state [Positive 
Affect and Negative Affect Scale, PANAS (Watson etal, 1988)], state 
anxiety level [Spielberger State and Trait Anxiety Inventory, STAI- 
State (Spielberger et al., 1983)], self-reported stress level [Perceived 
Stress Scale, PSS (Cohen et al, 1983)], trait rumination [Rumi- 
native Response Scale, RRS (Nolen-Hoeksema, 1991)], and trait 
mindfulness [Five Facet Mindfulness Questionnaire, FFMQ (Baer 
et al, 2006)]. Early-life stress was assessed using the Childhood 
Trauma Questionnaire (CTQ) (Bernstein and Fink, 1998), which 
contains five subscales: three assessing abuse (emotional, physical, 
and sexual) and two assessing neglect (emotional and physical). 
Early-life stress experiences were defined as either having at least 
three out the five subscales with a score >6 or meeting any of the 
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following cut-off scores (Bernstein and Fink, 1998): emotion abuse 
>9; emotion neglect >10; physical abuse >10; physical neglect 
>10; or sexual abuse >8. 

The imaging session was composed of an anatomical scan, a 
resting state scan, and four pairs of stress (or mindful breathing) 
task and an emotional go/no-go (EGNG) task runs. Using this 
design, we previously investigated the effect of stress and mindful- 
ness inductions on successful inhibition to negative vs. to neutral 
pictures during the EGNG task (Paul et al., 2013). In this report, we 
focus on the blood perfusion changes during while subjects per- 
formed the stress task and the mindful breathing task. Changes in 
heart rate, respiration rate, and Cortisol level were measured dur- 
ing the stress and mindful breathing tasks to objectively evaluate 
stress level. After completing each stress or mindful breathing task 
run, self-ratings of stress were also obtained immediately using a 
likert scale (ranged from 0 to 4, from not stressful at all to very 
stressful). Salivary Cortisol levels were measured at the beginning, 
middle, and end of the scan session. All imaging scans were com- 
pleted in the late afternoon to obtain low and stable Cortisol levels, 
when individuals are more responsive to stimulation (Jansen et al., 
1998) during these hours. Caffeine, smoking, and exercise were 
prohibited 2 h prior to scanning. 

EXPERIMENTAL DESIGN 

Stress induction task 

We used a mental arithmetic paradigm (Soufer et al., 1998; Wang 
et al., 2005) to induce stress similar to the Trier Social Stress Test 
(Kirschbaum et al., 1993). Participants were given a four-digit 
starting number and a two-digit integer (presented for 5 s) at the 
beginning of a run. Participants had to subtract a two-digit integer 
from the starting number and subtract continuously. The subtrac- 
tion was temporarily paused every 45 s when a fixation cross was 
presented, and continued when the fixation cross disappeared. Par- 
ticipants reported the final subtraction value at the completion of 
a run. Each run started with a different number and participants 
subtracted a different integer from the starting number during 
each run. Each run lasted 5 min. At the end of each run, subjects 
rated their stress level using a 1^ analog bar (with 1 being not 
stressful at all and 4 being very stressful) (Figure 1). 



Mindful breathing task 

During the pre-scan session of the mindful breathing task day, par- 
ticipants were trained to: (1) focus their attention on the bodily 
sensations of breathing and count breaths from 1 to 10, (2) notice 
their mind wander and return to counting when mind wander- 
ing happens, and (3) simply return attention to breathing without 
getting frustrated when their mind wanders. These instructions 
mirror a commonly used mindfulness meditation practice (Hanh, 
1976). Participants were instructed to follow these instructions 
during the mindful breathing task. For both stress and mindful 
tasks, participants paused when a fixation cross was displayed on 
the screen (Figure 1). 

We validated the stress and mindful breathing task in our previ- 
ous study. Subjects rated the stress induction task as more stressful, 
had higher Cortisol level, faster heat beat, and respiratory rate dur- 
ing the stress task across all runs than during the mindful breathing 
task (Paul et al., 2013). In this study, we investigated whether 
modulated moderated the effect of repeated mental stress. 

BIOCHEMICAL AND PHYSIOLOGICAL MEASURES 

To validate the stress and mindful tasks, peripheral biochemical 
and physical measures were collected including salivary Cortisol 
level, heart rate, and respiration rate. Participants were given a 
Salivette (Sarstedt AG & Co., Germany) to collect salivary right 
before, in the middle, and right after the whole scan session. 
Salivary Cortisol levels were assessed with solid-phase Coat-A- 
Count 125 I radioimmunoassays for Cortisol (TKCO) provided by 
Siemens Healthcare Diagnostics (Los Angeles, CA, USA) following 
our previous work (Schultheiss and Stanton, 2009; Stanton et al, 
2009). Heart rate and respiration rate were continuously mon- 
itored during scanning using a pulse oximeter and a chest belt, 
respectively (Biopac Systems, Goleta, CA, USA). 

IMAGE ACQUISITION AND ANALYSIS 

A 3.0 T GE MR750 scanner was used to acquire images at the 
Duke-UNC Brain Imaging and Analysis Center. After acquiring a 
Tl-weighted SPGR anatomical image (matrix= 256 x 256 x 180, 
1mm 3 ) axially parallel to the anterior and posterior commis- 
sures (AC-PC) line, we acquired perfusion images during a 5- 
min resting state and stress/mindful breathing tasks using an 



4p 



15~20s 



Subtract 
(orcount 



Subtract 
(or count] 



Stress 
rating 



Go 

► ► ► ► 



Nogo 



Subtract 
(or count 



stress (or relaxation) run- 
5min (ASL) 



Scr 



Neu 



_1_ 



t ► ► 

in 



3cr=Scramble: Neu=Neutral; Neg=Negative 

EGNG task run 

5min (BOLD) j 



One stress (or relaxation) - EGNG task pair. There were 4 pairs, i.e, repeated four times 

FIGURE 1 | Overview of the study design. There were four pairs of the stressful-EGNG tasks (also the same for the mindful-EGNG tasks). In this report, we 
focused our analysis on the stress (and mindful) task period and compared the repetition effect of stressful (and mindful) condition on perfusion changes over 
time (Run 1 vs. Run 4). 
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arterial spin labeling (ASL) sequence. A modified FAIR with spi- 
ral out sequence using a spatially selective inversion pulse and 
the QUIPSSII saturation pulses sequence (Wong et al, 1998) was 
used to allow quantitative determination of perfusion. A delay 
of 1 s was inserted between the end of the labeling pulse and 
image acquisition to reduce transit artifact. Acquisition parame- 
ters were: repetition time (TR)=4s, echo time (TE) = 3.2 ms, 
Til = 600 ms, TI2= 1600 ms, and flip angle = 90 °. Thirty-four 
slices (64 x 64 x 34 matrix; ~3.5mm 3 voxel size) were acquired 
from inferior to superior in a sequential order. BOLD images were 
acquired during the emotional go/no-go task. 

The analysis was carried out using FEAT-perfusion fMRI analy- 
sis (MRI Expert Analysis Tool Version 5.98), part of the FSL analy- 
sis package (FMRIB Software Library) 1 . We conducted full per- 
fusion signal modeling 2 , where three explanatory variables (EVs) 
were modeled: the BOLD signal, the (constant-height) tag-control 
difference, and the activation component of the tag-control signal 
(formed by multiplying the first two EVs together). Thus, the acti- 
vation component of the tag-control signal is our study interest. 
The following standard preprocessing steps were taken: removal of 
non-brain signal outside the head using the Brain Extraction Tool 
(BET), slice-time correction, coregistration, motion correction, 
normalization, spatial smoothing (5 mm FWHM), and high-pass 
filtering (1/60 Hz). 

Given that the primary goal of the study is to investigate the 
repetition effect of stress on perfusion changes over time and the 
modification effect of early-life stress, we focused our data analysis 
on between-run differences within subjects (i.e., the first run vs. the 
last run). Calculating the perfusion change from Run 1 to Run 4 is 
the way we modeled change over time, and in particular, response 
patterns to repeated stress. An increase in perfusion from Run 1 to 
Run 4 (i.e., Run 4-Run 1) demonstrates either neural sensitization 
in regions related to affective processing or more effective recruit- 
ment of resources in a brain region related to cognitive control, and 
a decrease from Run 1 to Run 4 demonstrates neural habituation 
or less effective recruitment of resources in a brain region related 
to cognitive control. The results were also confirmed by using the 
early two runs (Run 1 and Run 2) vs. the later two runs (Run 3 
and Run 4). The between-run differences were computed using a 
fixed effect model for each subject under both stressful task con- 
dition and mindful breathing task condition. For the group level 
analysis, we first examined significant perfusion between-run dif- 
ference in the group during the stress and mindful tasks separately 
using random effect model (FLAME 1) to investigate which regions 
demonstrated a habituation effect over the course of the stress task 
and the mindfulness task across all participants. Then, to under- 
stand whether the habituation effect was stress or mindfulness 
specific, we examined differences in neural adaptation patterns 
between the two tasks. We compared the difference in perfusion of 
Run 1 vs. Run 4 between the stressful and mindful task condition 
using paired (--test with a random effects analysis. 

To investigate the modification effect of early-life stress on 
habituation in response to repeated stress over time, we then 
conducted regression analyses to examine the association between 



1 www.fmrib.ox.ac.uk/fsl 

2 http://fsl.fmrib.ox.ac.uk/fsl/fsl4.0/feat5/perfusion.html 



perfusion changes from Run 1 to Run 4 with the measurements of 
early-life stress. Given the theory that early-life stress is associated 
with depression in adulthood, we first investigated the relationship 
between the total CTQ score and CTQ subcomponent scores with 
trait rumination (RRS) and trait mindfulness (FFMQ) as well as 
each facet of FFMQ using simple linear (total score of measure- 
ment) and multiple linear regression analysis (subcomponents). 
The total or subscores of CTQ which showed significant correla- 
tion (p < 0.05) with RRS or FFMQ were used to do the whole-brain 
voxel-wise analysis with perfusion changes from Run 1 to Run 4. 
For all imaging-related analyses, significance was determined using 
a voxel significance level of z > 2.3, with a whole-brain-corrected 
cluster significance threshold of p < 0.05. 

Finally, to further understand the perfusion pattern in response 
to the repeated stressor across time related to depression vul- 
nerability and resilience, we also subsequently identified indi- 
viduals with early-life stress who possibly were vulnerable or 
resilient to depression using the following criteria: (1) vul- 
nerable (n = 5), having early-life stress experiences, RRS > 40, 
FFMQ < 150, and (2) Resilient (n = 5), having early-life stress 
experiences, RRS < 40, FFMQ > 150. Individuals who did not 
have early-life stress experience were defined as Neutral (n = 6). 
The psychological measures for each of the groups were shown 
in Table 1. There were three participants who had missing RRS 
and FFMQ data. Because of the small sample size of the Vul- 
nerable, Resilient, and Neutral groups, we investigated the perfu- 
sion changes using Region-of-Interest (ROI) analysis. The ROIs 
were identified by extracting the significant clusters from the 
regression analyses in the third level analysis. The mean sig- 
nal strength with each ROI for each run in each subject was 
calculated. 

RESULTS 

BEHAVIORAL RESULTS 

The negative association of early-life stress with trait mindfulness 

There was not a significant correlation between RRS and either the 
CTQ total score or the score of subcomponents. However, multiple 



Table 1 | The Neuropsychological measures for the three groups, 
group resilient to depression, the group vulnerable to depression, and 
the group neutral to depression that were defined based on the RRS 
and FFMQ scores. 



No. of subjects 


Resilient 


Vulnerable 


Neutral 




5 


5 


6 


RRS 


26.4 (3.2) 


43.6 (3.0) 


28.5 (4.0) 


FFMQ 


152.2 (3.5) 


125.2 (11.9) 


115.3 (16.5) 


Total CTQ 


28.6 (3.3) 


32.2 (3.3) 


29.8 (6.5) 


Emotional abuse 


5.2 (0.4) 


8 (1.4) 


7 (2.8) 


Physical abuse 


5.8 (1.3) 


5.6 (0.9) 


5.5 (0.8) 


Sexual abuse 


5(0) 


5.4 (0.9) 


5(0) 


Emotional neglect 


7 (1.6) 


7.8 (2.9) 


6.8 (2.6) 


Physical neglect 


5.6 (0.9) 


5.4 (0.5) 


5.5 (0.8) 


Italics are used when referring to the subscales of the CTQ (Childhood Trauma 



Questionnaire). 
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regression analysis showed that, among the CTQ subcomponents, 
emotional abuse was inversely associated with non-reactivity, one 
of the five facets of the mindfulness on the FFMQ {t = 2.23, 
p = 0.05), meaning that individuals with more severe emotional 
abuse during childhood had lower trait non-reactivity. Trait non- 
reactivity is the ability to step back from one's experiences without 
becoming overly engrossed in them (Baer et al., 2006), which we 
previously found having protective effect from depression vulner- 
ability (Paul et al., 2013). Therefore, emotional abuse seemed to 
be an important risk factor for depression vulnerability, and the 
score of emotion abuse was used to our following neuroimaging 
data analysis. 



Stressful Task Mindful Breathing Task 







n 





FIGURE 2 | The main effect of task repetition (Run 1-Run 4) on 
perfusion during the stress (left) and mindful breathing (right) tasks 
(voxel-wise analysis, Z > 2.3, p < 0.05 cluster correction). HC, 

hippocampus cortex. 



NEUROIMAGING RESULTS 

The habituation to the repeated stressor and repeated mindful task 

First, we investigated regions that showed a habituation effect, i.e., 
regions that showed perfusion decrease from Run 1 to Run 4 (Run 
1-Run 4) under the stress and mindful breathing conditions sepa- 
rately across all participants. The voxel-wise analyses revealed that 
when a stressor was presented repeatedly (Figure 2; Table 2), the 
perfusion was significantly decreased in the left middle and supe- 
rior temporal cortex from the first to the last run (i.e., Run 1-Run 
4) suggesting there was a habituation effect in these regions under 
the stressful condition. During the repeated mindful breathing 
task (Figure 2; Table 2), the perfusion was significantly decreased 
in the left amygdala, hippocampus, fusiform gyrus, and bilateral 
thalamus, therefore, there was a habituation effect in these regions 
during the mindfulness induction. In a direct comparison of the 
stress task and mindful breathing task, however, we did not find 
significant differences in perfusion comparing the Run 1-Run 4 
contrast between the 2 days. 

The association of early-life stress with less habituation to the 
repeated stressor 

Given that only the emotional abuse subscore of the child trauma 
score was significantly correlated with mindfulness (negatively), 
we used the emotional abuse subscore to compute the regression 
analyses with the perfusion change from Run 1 to Run 4. The 
whole-brain voxel-wise analysis revealed that the emotional abuse 
subscore was significantly correlated with a reduction in perfu- 
sion level from Run 1 to Run 4 in the right orbitofrontal cortex 
and the vmPFC (Figure 3; Table 3) during the stress task. Taking 
the significant cluster of vmPFC as a ROI, we confirmed that the 
significant correlation between the reduction of perfusion from 
Run 1 to Run 4 and emotional abuse in childhood (ru = 0.59, 
p = 0.016, for double confirm purpose) was not due to outliers 
(see the lower plot in Figure 3). The post hoc ROI analysis showed 
that the perfusion reduction from Run 1 to Run 4 during the 



Table 2 | Regions that revealed habituation effect (Run 1-Run 4) of the stressful task and the mindful breathing task (those in italic are 
subclusters). 



Region 



Hemisphere 



Cluster size 



Peak voxel (x) 



Peak voxel (y) 



Peak voxel (z) 



Maximum Z-value 



13957 



HABITUATION EFFECT OFTHE STRESSFULTASK 

Cerebellum cms I 
amygdala 
hippocampus 
temporal pole 



R 
L 
L 
L 
R 
L 
L 
R 

R 1801 
HABITUATION EFFECT OFTHE MINDFUL BREATHING TASK 

Midbrain L 1862 

midbrain Ft 
pons Ft 
thalamus L 
hippocampus L 



cerebellum vermis 
fusiform gyrus 

Lateral parietal cortex 



28 
-28 
-32 
-42 

40 
-6 
-28 

36 

32 

-6 

6 

2 
-20 
-40 



-80 
-6 
-26 
-12 

0 
-49 
-56 
-36 
-70 

-16 
-10 
-20 
-20 
-20 



-22 
-14 
-11 
-30 
-36 
-23 
-6 
19 
54 

-16 
-12 
-24 
-10 
-16 



4.36 
3.74 
3.01 
4 

3.92 
3.28 
3.52 
2.87 
3.47 

4.7 

3.7 

4.25 

3.71 

3.73 
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stressful condition was also significantly correlated with the CTQ 
total score (ru = 0.56, p = 0.024, Figure 3). A similar association 
was not found in response to the repeated mindful breathing task, 
indicating that reported early-life stress was uniquely associated 
with a habituation response in the vmPFC perfusion under acute 
stress conditions in the laboratory. 

Habituation effect in depression-vulnerable vs. resilient 
individuals - subsequent ROI analyses 

We further performed an exploratory investigation on the perfu- 
sion pattern in the vmPFC and in the amygdala in response to 
the repeated stressor in depression-vulnerable and resilient indi- 
viduals. The vmPFC ROI was defined using voxels that revealed 
significant correlation between emotional abuse and the reduc- 
tion of Run 1-Run 4 perfusion, and the amygdala ROI was 
defined using voxels showing significant habituation effect during 
the mindful breathing task. Only ROI analyses were conducted 
due to the small number of subjects who met the vulnerable 
or resilient criteria. We focused on the vmPFC and amygdala 
regions because of the well-known association of these two regions 
with early trauma (Taylor et al., 2006; Hart and Rubia, 2012) 
and the inverse relationship between the two regions in emotion 
regulation. 



As the bar graph shown in Figure 4, ROI analyses on the vmPFC 
and amygdala revealed an inverse relationship. For individuals vul- 
nerable to depression, the perfusion level in the amygdala increased 
from Run 2 to Run 4 relative to Run 1, whereas the perfusion level 
in the vmPFC was decreased from Run 3 to Run 4 relative to Run 
1. Interestingly, depression-resilient individuals showed a linear 
decrease in perfusion level from Run 1 to Run 4 in the amygdala, 
whereas the perfusion level of the vmPFC did not change signifi- 
cantly across the four runs. By contrast, individuals who did not 
experience childhood trauma (N in Figure 4) revealed a sustained 
activity perfusion level in the vmPFC when they were repeatedly 
exposure to stress and sustained activation in the amygdala. 

DISCUSSION 

The current study investigated the impact of early-life stress on 
blood flow over time during a repeated stress induction task and 
during a repeated mindful breathing task. We found a inverse 
correlation between emotional abuse and non-reactivity (a sub- 
component of mindfulness) which is consistent with previous 
finding that childhood emotional abuse and neglect proved more 
predictive of adult depression than childhood sexual or physical 
abuse (Powers et al, 2009). Our voxel-wise regression analysis 
revealed a clear association between childhood emotional abuse 
and reduced orbitofrontal and vmPFC perfusion over time when 
participants were exposed to stress repeatedly. The orbitofrontal 
and vmPFC sends projections to the amygdala, and is known 
to regulate amygdala activity (Price, 2005; Urry et al., 2006). 
Therefore, we speculate that a lack of sustained orbitofrontal and 
vmPFC activity, coupled with amygdala hyperactivity, in response 
to a repeated stressor, might be an important neural signature of 
depression vulnerability. 

Our interpretation was supported by our subsequent ROI 
analyses which confirmed that individuals who experienced early- 
life stress and were high in trait rumination had reduced vmPFC 
activity in the later compared to earlier stressful trials, individuals 
who experienced early-life stress but were high in trait mindful- 
ness had sustained vmPFC activity, and individuals who did not 
have early- life stress had an increased vmPFC perfusion over time. 
The inverse relationship between the vmPFC and the amygdala 
has been consistently reported in the literature (Phelps et al, 2004; 
Harenski and Hamann, 2006; Ohira et al, 2006; Urry et al, 2006). 
Our results suggest that healthy adults exert effortful control over 
emotion via increasing the vmPFC activity when they are repeat- 
edly exposed to a stressor. For individuals who develop stress 
resilience, they may regulate amygdala activity more efficiently 
than healthy adults without this trait resilience by sustaining the 
vmPFC activity without further increasing the effortful vmPFC 
activity. Whereas when individuals are vulnerable to stress, the 
vmPFC has a limited capacity for cognitive control to regulate 
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FIGURE 3 | Left, voxel-wise analysis on the correlation between 
emotional abuse score and the perfusion changes from Run 1 to Run 4 
during repeated stress (voxel-wise analysis, z > 2.3, p < 0.05 cluster 
correction). Right, correlation plots based on ROI analyses confirming the 
significant correlation between the perfusion changes from Run 1 to Run 4 
during repeated stress and emotional abuse score (lower plot) and 
ChildhoodTrauma Questionnaire (CTQ) Total Score (upper plot). 



Table 3 | Regions that revealed significant correlation between perfusion reduction during repeated stress and emotional abuse (those in italic 
are subclusters). 

Region Hemisphere Cluster size Peak voxel (x) Peak voxel (v) Peak voxel (z) Maximum Z-value 

Orbital frontal cortex R 1909 36 54 -10 4.68 

Orbital frontal cortex R 34 60 -4 3.84 

Ventromedial prefrontal cortex M -1 45 -6 3.25 
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FIGURE 4 | The perfusion changing pattern in the vmPFC and amygdala 
to repeated stress across runs among individuals resilient to depression 
(/?), vulnerable to depression (V), and neutral [N, not resilient and not 
vulnerable). Red = Run 1 ; pink = Run 2; silver = Run 3; blue= Run 4. The 
perfusion level in the amygdala was increased from Run 2 to Run 4 relative to 
Run 1, whereas the perfusion level in the vmPFC was decreased from Run 3 



emotional responding, and could not effectively regulate amygdala 
activity in the later runs. Therefore, our study suggests that a lack of 
sustained vmPFC activity when repeatedly exposed to stress might 
be an important neuroimaging marker that distinguishes stress 
vulnerability vs. resilience in individuals with early-life stress. 

However, the interpretation is preliminary due to the small 
sample size for the stress vulnerable and resilient groups. Another 
caveat of the study is the fact that we only recruited males in the 
study. The literature on sex differences in stress response is mixed 
(Kudielka and Kirschbaum, 2005; Wang et al., 2007; Ossewaarde 
et al., 2010). According to Ossewaarde et al. (2010), women have 
different responses to stress and facial expressions during different 
menstrual phases. Their results also indicated general stress sen- 
sitivity during the late luteal phase, as indicated by stress-induced 
heart rate changes and reported negative affect. Indeed, a study by 
Kirschbaum et al. ( 1993 ) found that the salivary Cortisol response 
of women in the luteal phase was similar to men, but women 
taking oral contraceptives or in the follicular phase demonstrated 
lower free Cortisol responses. Wang et al. (2007) further showed 
that stress in men was associated with cerebral blood flow (CBF) 
increase in the right prefrontal cortex and CBF reduction in the 
left orbitofrontal cortex, whereas stress in women activated the 
ventral striatum, putamen, insula, and cingulate cortex. There, 
future studies with larger sample sizes of both males and females 
are needed to verify our findings. Studies include women, and 
either collect menstrual-cycle phase data or only examine women 
in the luteal phases are recommended, to minimize the effect of 
menstrual-cycle phase on stress reactivity and further clarify the 
literature on sex-related differences in stress reactivity. In addi- 
tion, future studies using causal modeling are necessary to further 
understand whether a lack of regulatory control by the vmPFC 
over the amygdala is the most important regulatory pathway in 
depression vulnerability. 

Of note, although the hippocampus is often regarded as the 
most vulnerable region to stress and has consistently been found 
damaged in patients with PTSD (Astur et al, 2006; Werner et al, 
2009; Dickie et al, 2011; Thomaes et al, 2013), we did not find 



to Run 4 relative to Run 1 for individuals vulnerable to depression. On the 
other hand, depression-resilient individuals showed a linear decrease in 
perfusion level from Run 1 to Run 4, whereas the perfusion level of the 
vmPFC did not change significantly across the four runs. "Indicates significant 
difference between the first two runs with the last two runs in ROI analyses 
(paired f-tests, p<0.05). 



an association of perfusion change over time in the hippocampus 
with early-life stress. Instead, we found reduced perfusion in the 
hippocampus when repeatedly practicing mindful task suggesting 
a task learning effect. The reason that we did not find a significant 
association between the perfusion change in the hippocampus and 
early-life stress might be due to the fact that some individuals with 
early-life stress in our sample have developed resilience to stress. 
Our subsequent analysis with the amygdala demonstrated the pat- 
tern we would have expected with the hippocampus, and future 
studies in a larger sample size would clarify this hypothesis about 
stress and the hippocampus. 

We also found decreased blood flow over time in the amyg- 
dala, fusiform gyrus, and bilateral thalamus specifically during the 
mindful breathing task. The reduced blood perfusion from Run 
1 to Run 4 under the mindfulness condition fits well with previ- 
ous knowledge in that repeated mindful breathing can result in 
habitation and less excitation over time in the regions related to 
emotional salience (amygdala, thalamus) and imagery memory 
(hippocampus and fusiform gyrus) (Byrne et al, 2007; Bird et al, 
2012). Nevertheless, we did not find any significant difference in 
the perfusion pattern in the vmPFC or any other regions when 
comparing perfusion patterns over time between stress and mind- 
ful conditions. One explanation for this result might be that the 
stress task did not evoke a sufficiently strong stress response to 
see a statistical difference in our small sample size. Therefore, only 
individuals who had experienced early-life stress revealed stronger 
change in neural response to the repetition of such a mild stres- 
sor, suggesting that early-life stress might have a sensitized effect 
in these individuals. In this sense, a mild stressor might be more 
useful in detecting stress vulnerability. On one hand, our mental 
arithmetic task only involved mental effort but not negative feed- 
back or psychological pressure, which are two important elements 
of mental stress (Dickerson and Kemeny, 2004). It will be infor- 
mative in future studies to compare neural responses to different 
elements and different types of stressors. 

In summary, we found an association between early-life stress 
and reduced perfusion level in the vmPFC in response to repetition 
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of a mild stressor. The effect was particularly strong in indi- 
viduals who had early-life stress and also developed trait rumi- 
nation, but not in those who had strong trait mindfulness. 
Our results open an opportunity of assessing stress vulnerabil- 
ity and resilience by continuously monitoring the vmPFC and 
amygdala perfusion or activation change over time. Functional 
MRI typically has large variation across subjects and difficult to 
establish a norm range, which has limited its application clini- 
cally for individualized diagnosis. The dynamic changing pattern 
in brain perfusion over time found in this study is important 
because it relies on the changing pattern within subject across 
time rather than a norm value. If confirmed by future studies 
using larger samples of males and females, changing patterns of 
perfusion to repeated exposures of stress over time, could be a 



novel way to identify depression vulnerability at an individual 
level. 
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